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Abstract 

We consider the most general renormalizable chiral Yukawa model with SU(3) co i OI 
replaced by SU(N C ), SU(2) L replaced by SU(N W ) and U(1) Y replaced by ^(l)^" 1 in 

the limit N c — > oo, N w — ► oo with the ratio p = ^f- ^ 0, oo held fixed. Since for N w > 3 

only one renormalizable Yukawa coupling per family exists and there is no mixing between 
families the limit is appropriate for the description of the effects of a heavy top quark when 
all the other fermions are taken to be massless. A rough estimate of the triviality bound 
on the Yukawa coupling is equivalent to m t < 1 TeV. 
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In the minimal standard model at energies above 100 GeV only the scalar selfcoupling 
A is completely unknown and might be strong. Skeptics 1 might argue that a strong Yukawa 
coupling y has not been yet ruled out completely and this note addresses the problem how 
to do calculations if they are right [1]. 

Excepting the A and y all other couplings are weak and can be treated perturbatively. 
The problem then is to set up a computation first in the case where all the weak couplings 
are set to zero. We have some requirements: Because of "triviality" and "vacuum stability" 
issues one ought to be able to keep an explicit cutoff in the calculations. The cutoff has 
to be defined in a reasonable manner, meaning that an expansion in inverse powers of the 
cutoff can be set up and has the normal structure one would expect of terms induced by a 
higher embedding theory (for example, a sharp momentum cutoff is not acceptable). We 
do not want to use the lattice since many of the non-perturbative effects in lattice models 
have little to do with continuum physics. 2 Intuitively we feel that a sensible approximation 
should keep more or less a democratic view of fermionic and scalar excitations. 

We propose that a good way to achieve the above is to consider a model with a 
single Higgs multiplet where SU(3) co \or is replaced by SU(N C ), SU(2)l is replaced by 
SU(N W ) and U(l)y is replaced by ^(l)^ -1 in the limit N c — > oo, N w — ► oo with the 



ratio P = \j yn~ ^ ^' 00 ne ^ fi xe d- Just letting N C (N W ) go to infinity with N W (N C ) fixed 
suppresses bosons (righthanded fermions) and these limits may distort effects like vacuum 
destabilization [2] or bag formation [3]. In any case more diagrams get summed up when 
both N c and N w go to infinity and all diagrams you would get when only N c or N w go to 
infinity are included. 

We introduce the model in Euclidean space; with a finite cutoff the model is non- 
perturbatively well defined but the continuation to Minkowski space would introduce uni- 
tarity violating effects at energies of the order of the cutoff. This drawback is compensated 
by preserving translational and full rotational invariance with propagators regulated in a 
Pauli-Villars manner. The action is given by S: 



#? (r) M? (r) + xt )s h x (-d 2 W-x- {r)s + ^M-aW - - 4]^*<« 2 

+4f £0- (r) x- r,) >'4S + xt )s ^ {r,) ^KV] 



1 The assumption that the top Yukawa coupling is relatively weak is consistent with 
present knowledge indicating that m t <~ 200 GeV . 

2 The lattice introduces in a smooth manner a finite bottom to the Dirac sea and this is 
particularly disturbing to chiral fermions. One could however use our technique to study 
non-perturbative effects in lattice models and hopefully sort out the situation there. 
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where 

N =VN W N C , 0+ = a M <9 M , 0. = 

ai=<Jl, (72 = &2, 0"3 = 0"3, 0"4 = — 0"4 = — « (2) 

z=l,---,iV w a=l,---,iV c r=l,---,nf s = l,---.N w 

The £Ti 2,3 are Pauli matrices. The functions h x and h<j> can be chosen in many ways and 
introduce the regularization. Ultraviolet divergences can be cured by just regularizing the 
X and <p propagators because there are no loops consisting of only ip propagators. When 
in need of an explicit form we shall use h x (p 2 ) = h ( / ) (p 2 ) = 1 + (p 2 /A 2 ) n with a sufficiently 
large n. 

The symmetry properties are as follows: The space-time group is SU(2) x SU (2) and 
a spinor with subscript + transforms as (|,0) while one with subscript — transforms as 
(0, 2)- ip, Xi X are two component Grassmann spinors and complex conjugation does not 
act on them. Under SU(N C ) x SU(N W ) the representations are t(j- ~ (N C ,N W ), ip+ ~ 
(N C ,N W ), x+ ~ (N c ,l), X- ~ (^,1), <P* ~ (l,N w ), (f> ~ (1,JV W ). This representation 
content is independent of the indices r and s whenever they appear; the index r runs 
over families and s over the number of fermionic weak singlets (and color mutiplets) per 
family. The representations under the remaining t/(l)'s are chosen subject to anomaly 
considerations when gauging is contemplated and eliminate all but one of the matrices 
A^ s \ The non-vanishing matrix is chosen to be and the remaining x fields with 
index s > 2 decouple at zero gauge couplings. When N w = 2 an additional set of Yukawa 
couplings is allowed by SU(N C ) x SU(N W ) and also by the t/(l)'s. These couplings do 

not generalize to N w > 3 and will henceforth be ignored. The matrix A^J, can be made 
diagonal with positive entries by a bi-unitary transformation, decoupling the families. In 
conclusion one can ignore the r and s indices, which is just as well because we want to 
keep only one potentially heavy "top" color multiplet. The model then becomes [1]: 



-S = [ V!_^ + V; a + X a +h x (-d 2 )0-X a - + <TM-<9 2 )dV - mji^V* - -^-(0**0* ) 2 
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When iV — > 00 with p = \J^f~ held fixed the dominating diagrams are planar. Renor- 

malizability allows us to get away with not having to introduce four Fermion interactions 
and with them an insoluble planar diagram problem, at no cost to generality. 3 Thus the 
model is manageable being similar diagrammatically to two dimensional QCD with matter 
at large N c . The large iV limit can be found by looking at the structure of planar Feynman 
graphs: essentially one sums over all "cactuses of bubbles" and all "rainbows" . For future 
work it is somewhat preferable to use functional integral manipulations to the same end. 

In this note we wish to show feasibility within a simplest nontrivial example. We set 
our modest goal to compute a /9-function associated with the Yukawa coupling constant 
defined in some simple manner and ignoring cutoff effects. This can be done with relative 



Except for possible cutoff effects. 
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ease because even at infinite order in the couplings the flow of this Yukawa coupling is 
unaffected by Ao, similarly to the finite N form at one loop order. From this calculation a 
rough estimate for the "triviality" bound on the top mass can be obtained. 4 

To find (3 functions it suffices to solve the model at criticality when all masses vanish. 
We start in the symmetric phase and pick a mass independent renormalization scheme to 
make the zero mass limit smooth. 

We introduce bilocal fields, 

SaA*,v) = ^X a +a {x)x%{y), K(x,y) = ^WG/), (4) 

where a, (3 are space-time spinor indices, into the functional integral by writing representa- 
tions for the appropriate 5-functions with the help of two auxiliary bilocal fields A a ,p(x, y) 
and fj,(x,y). We can integrate out now all the original fields, including the zero mode of 
the scalar field because we assumed that we are in the symmetric phase. This makes the N 
and p dependence explicit and, when iV — > oo, the integral over S, A, K and fx is dominated 
by a saddle point which has to satisfy the following equations: 

K(x, y) ={[-M-d 2 )<9 2 + ml+ fA'^x 
S a ,p{x,y) ={[h x (-d 2 )0- +\]~ 1 } yl 3, xa 

\xA x > y) = ~ 9 2 p^{0+ 1 }xa, y pK(x, y) (5) 

' a, P 

We go to Fourier space and use translational and rotational invariance to write 

KAp) = ip^^Hp 2 ) Kp) = Kp 2 ) (6) 

Let us introduce Z x and wave function renormalization constants, and new functions 
/ = Z x [X(p 2 ) + h x (p 2 )] and g(p 2 ) = Z ( j } [p 2 h ( j ) (p 2 ) + m\ + fi{p 2 )]. Since at the planar 
level there are no vertex corrections to the Yukawa coupling 5 it makes sense to define a 

2 

renormalized Yukawa coupling by g\ = Z^Z x g\, (a = -jf^)- We also set A' = XqZ 2 but, 
unlike a, do not expect X' Q to stay finite when A — > oo. Z x and are defined in terms of 

the bare couplings by f(fi 2 ) = g ^ 2 - = 1- 

After performing the angular part of the momentum integrals one gets when the 
products in the saddle point equations for x, y dependent quantities are converted into 



4 In practice this bound might be somewhat uninteresting because it is likely to be 
higher than the vacuum stability bound for reasonable Higgs masses. 

5 This is also true at one loop order at finite N. 



4 



convolutions in momentum space one finds: 
f(p 2 ) =Z x h x (p 2 ) + V2ap 



p k 4 dk 2 

2 J p2 g(k 2 ) ' p 2 J g(k 2 ) 2p* J ~gj&) 



i f°° dk 2 i r k 2 dk 2 

+ 



g(p 2 ) =Z <l> p 2 h 4> {j?) + 



a 



p2 r°° dk 2 r p " dk 2 i r 

~2 1 P * WW) + Jo P Vio 



; k 2 dk 2 

1W) 



+ m 2 (7) 



X' r» k 2 dk 2 V2a f 



dk 2 
f(k 2 ) 



The massless case is obtained by adjusting TOq so that m 2 = and then it makes sense to 
define h(p 2 ) = 9 ? ' leading to a set of equations of a more symmetric appearance: 



h(p 2 ) =z^(p 2 ) + y^c[f]( P 2 ) 

f(p 2 ) =Z x h x (p 2 ) + V2ap£[h](p 2 ) 



i r 



dk< 



1 



dk< 



i r k 2 dk 2 



2 k 2 X(k 2 ) p 2 J X(k 2 ) 2p 4 J X(k 2 ) 



(8) 



f(p 2 ) =h(p 2 ) = 1 

Note the interesting bosonic-fermionic symmetry at p = 1. We now rescale the functions 
and the momentum variable by: 



f(up 2 ) = v2a(f){u) h(up 2 ) = ^faip{u) 



(9) 



and obtain: 



4>{u) 



1 4[mV)-m. 2 )]^ fUv r 1{1 - x)dx 



a 
1 



+ 



v Jo 



X 



2^/a 2^/a 



[h x (up 2 ) - h x (p 2 )] +p 



4>{xv) 
1 dv f 1 (1 — x)dx 
v J tfj(xv) 



(10) 



It is clear that when n is large enough h and / are dominated by the h functions in 
the ultraviolet and therefore the large cutoff behavior of the wave function renormalization 
constants is such that in the infinite cutoff limit at u and p 2 fixed Z^ and Z x simply 
disappear from the above equation leaving us with 



1 : 1 r 1 dv f 1 (1 - x)dx 
v 



4>{u) = 



l 



2v^ 



+ P - 



(j)(xv) 
1 dv f 1 (1 — x)dx 
^(xv) 
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On physical grounds it is obvious that we wish that V> and be positive. The equa- 
tions above show that if this is true in some interval (0, it*) then the functions will be 
monotonically decreasing there. However, the equation do not admit an asymptotic be- 
havior as u — > oo with both functions approaching non-negative limits so the positivity 
requirement must get violated somewhere in the ultraviolet. 6 Starting at the scale where 
the violation first occurs cutoff effects cannot be neglected any more. This "unphysical" 
scale is the usual Landau pole, this time appearing in a nonperturbative approximation. 

In the infrared there is no Landau pole problem and cutoff effects are indeed negligible. 
To investigate the behavior there we set x = — log (it) and tfj(x) = ip(u), 4>{x) = 4>{u) and 
derive: 

^0=7 (12) 

ip 

V-^ 3^ 2 d - d ( d l\( d 2 

dx 3 dx 2 dx dx \dx J \ dx 



These equations allow us to extract the infrared behavior of the solution of the integral 
equations. We ended up with only differential equations (rather than integral) reflecting 
that according to the Renormalization Group one needs only very limited information at 
a given scale in order to derive the behavior at a scale close by; therefore a differential 
equation must show up eventually, its order less the number of asymptotic conditions 
determining the amount of information at a given scale that is needed to determine the 
behavior at the next scale. Our choice of regularization was made so that even at finite 
cutoff <p and i/j obey a set of purely differential equations of a structure similar to (12). 

Since the fixed point governing the infrared behavior is the free field fixed point the 
form of the solutions in the infrared simply embodies the two anomalous dimensions as- 
sociated with the x an d 4> fields when expanded in a around a = 0. These anomalous 
dimensions as a function of the coupling also determine the (3 function and from the asymp- 
totic series of i/j(x) and 4>{x) at x — > oo we can get the contributions to the (3(a) function 
ordered in the number of loops. 

The differential equations lead to: 

^~^ 6 [1+C^ +€'- + •••] 

X X 

t-fox^l + d^ + d '- + ...] 

XX 



(13) 



6 A zero of ip or at a positive u corresponds to a pole in a two-point function at an 
Euclidean momentum; the "particle" associated with this pole would be tachyonic. 
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where 

a=hp =TT? 



d = cp2 = -j^w (14) 



d> , 3(1 - p 2 ) p + p- 1 

= 2(1 = — ^— 

Only the product V'o^o an d a particular linear combination of d! and c' get determined 
because the equations are invariant under a field rescaling ip — * A -1 -!/*, — * reflecting 
a change in the finite parts of the wave function renormalization constants and a shift 
x — > a; + xo reflecting a change in p 2 . The field-rescaling invariance disappears in the 
product r(u) = v / 2 , 0('u)</>('u) and this combination is indeed special because it satisfies a 
Renormalization Group equation in which only the /3-function appears (the sum of the 
two anomalous dimensions rather than each individual one): 



r( e 2t ,a)=0 (15) 



This equation is exact in our limit. From it we derive another exact relation between 
r and (3: 



2u _ or„/„ . 2m2 dr ( U ) 



f3(a(up )) = — 2[a(up z 



d(logu) 



(16) 



The argument of a in the equation above can of course be ignored. 

Using now the asymptotic expansion of the exact solution in the infrared we get: 

f3(a)=a 2 [V2(p + p- 1 )-3a + 0(a 2 )} (17) 

The one loop result when expressed in a conventional variable, y(t) (in the standard model 
the top mass at tree level is m t = y(~ 0) 246 [GeV]), is, with the true value p = 

= ^2 which is close to the finite iV result, — 1^. The difference comes from the 

at tt z ' at 8ir z 

wave function renormalization of the lefthanded ip field which is suppressed at leading order 
in jj. For arbitrary iV c and N w we would have gotten at one loop order = ( 2N ^+^+ 1 )y 
and we see that the relative magnitude of the 1/N correction is much smaller than in the 
scalar case [4]. 

To get the full /3-function at infinite iV we solve for i/j and <p numerically by iterating 
an integral form of the equations based on (7) but with the normalization conditions at p? 
incorporated; this avoids the appearance of double integrals like in (11-12). The procedure 
is probably safer than using the third order differential equations and trying to enforce 
boundary conditions at x — > oo that eliminate the exponentially growing components e x 
and e 2x (see the last line in equation (12)). The result for (3 is shown in the figure and 
one again sees explicit evidence for the nonperturbative existence of a Landau pole. It 
is educational to compare the Landau pole energy obtained from the exact result to the 
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one loop estimate. For example, with a(/j 2 ) = 1.365 we get AL andau = 2/j 2 while at one 
loop we get AL andau = 1.565/U 2 . The value of ct(/i 2 ) we chose is very high because cutoff 
effects ought to be substantial when \i is so close to ALandau- A top quark corresponding 
to such a strong coupling would have a mass of rrtt = ^\p^ 246 [GeV] ~ 1043 GeV. 
We see that "triviality bounds" on the top will likely come out close to the perturbative 
unitarity bounds of Chanowitz et. al. [5] a property of our approximation that we view 
as an improvement on the work of Einhorn and Goldberg [6] who obtained 5600 GeV. 
Our estimate is close to that obtained from simple variants of large N W (N C fixed) or large 
N C (N W fixed) expansions [7]; this agreement holds essentially because a truncation of the 
set of planar diagrams to one loop turns out to be relatively acceptable numerically. From 
our non-perturbative result in the figure one sees that the one and two loop errors in j3 are 
about comparable for a ~ 0.4 (roughly of order 20 percent) and that using the one loop 
results all the way out to infinity is not very harmful, but the two loop result would be 
completely misleading if extrapolated beyond its region of validity, to a ~ 1. A top mass 
of 200 GeV would correspond to a m 0.05, well within the perturbative domain. 

We hope to work out a more comprehensive analysis of the N = oo limit in the future, 
starting with a calculation of a /3-function associated with Ao- 
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Figure Caption: The nonperturbatively determined /3-function for the Yukawa coupling 
at iV c , N w — > oo with = |. At tree level the top quark mass is given by m t = 

47r^/y| 246 [GeV]; /3(a) = ^ where t is the logarithm of the energy scale. 
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